Production and Detection of Randall-Sundrum 
Light Radions at Hadron Colliders 



Uma Mahanta^ and Anindya Datta*^ 
Mehta Research Institute, Chhatnag Road, 
Jhusi, Allahabad-211019, India 

e-mail; ^ mahanta@mri.ernet.in, anindya@mri.ernet.in 



Abstract 

In this paper we use the conformal anomaly in QCD to derive the coupling 
of light radion to gluons in the Randall-Sundrum model and use it to compute 
the radion production cross section at hadron colliders by gluon fusion. We find 
that the radion production cross section by gluon fusion at LHC would exceed 
that of the higgs boson by a factor that lies between 7 and 8 over most of the 
range. The decay modes of the radion are similar to that of the SM higgs boson. 
But the striking feature is the enhancement of radion to 2-photon and radion 
to 2-gluon branching ratio over the SM case. Utilising this, we then discuss the 
possible search strategies of such scalars at Tevatron and LHC. Using the 77 
decay mode one can explore/exclude radion mass upto 1 TeV. Even with the 
current collected data at the Tevatron, one can exclude radion mass upto 120 
GeV for {(j))= 1 TeV. 

PACS number (s): 



1 Couplings and Decay of the Radion to the SM 
particles 

Recently Goldberger and Wise showed that the modulus in the Randall-Sundrum 
[H scenario can be stabilized Q by introducing a scalar field in the bulk. The 
stabilized modulus turns out to be very light if its mass arises from a small bulk 
scalar mass. It was shown subsequently that the radion couples to the SM 
fields on the visible brane via the trace of the energy momentum tensor. In this 
article we shall study the production and decay of such a light radion at hadron 
colliders. 

The radion couples to the SM particles on the visible brane as given in M 
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via the relation, Lint — -^'rii4> where is the trace of the symmetrized and 
conserved energy momentum tensor for SM fields. At the tree level it is given 

by, 

= E + 2Mf Z^Z^ + 2ACW''Wf, + 2mlh\ (1) 

/ 

The fermion and gauge boson terms show the scale breaking effects due to 
electroweak symmetry breaking. Since the gluons and photons are massless the 
radion does not couple to these at the tree level. However the running of the 
gauge coupling in QCD and QED breaks the scale invariance and induces a 
trace anomaly [Q. The trace anomaly in QCD therefore generates the radion 
coupling to gluons which is given by 



where ^ - -(11 - ^Uf)^. G", = M - d,g- + gP'^g'^gt] is the 
gluon field strength tensor. For m?- < 4mj we have Uf — 5 dynamical quarks and 

hence SMA ~ _3_84Hs Qn the other hand for m^- > 4™? we have Uf — 6 and 

2gs 4:11 (f) ^ ■> 

hence Sis^l ~ —3.50^. The anomaly contribution is independent of fermion 
mass. Even if EW symmetry were exact and all fermions had remained massless 
the trace anomaly would still lead to the above coupling of to gluons. In the 
presence of EWSB, heavy quark loops give rise to a contribution to C^^g. In 
the infinite mass limit this contribution can be obtained by replacing v by (0) 
in the effective Lagrangian for hgg coupling [eqn. 10]. It can be shown that 
this contribution is smaller than the anomaly contribution written above. The 
running of QED coupling also introduces a conformal anomaly. This gives rise 
to the following coupling of the radion to the photons. 



{cf>) 2e 



Where F'^'^ has the usual meaning. We have used the following values of 
in our calculation. 
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2mt > > 2mw 

mi < 2mw (4) 



We want to mention that in the SM, in the heavy quark limit, the /177 
coupling is given by, 
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r^^^^^lM/ji^ F'^- (5) 

V 2e 



where, 

/3(e) 2a 



(6) 



2e Qtt 

Hence for m^, m^^ < 4 (where — > 77 is significant and the heavy top limit 

is vahd) we have = 57^. Although (j)ff and (j)VV iV = W, Z) couphngs 

are suppressed relative to hff and hVV couplings for {(j)) — 1 TeV, (pgg and 
(5(177 couplings are slightly enhanced relative to hgg and ft.77 couplings. 




Figure f: Branching ratio of Radion into different channels for {(j)) = 1 TeV 

Now we are in a position to present the radion decay branching ratios into 
different channels, fn fig. we present the relevant decay branching ratios. For 
the illustrative purpose we assumed {(p) = f TeV. The radion coupling to any 
SM field has the -t^ dependence. Thus if we choose a different value of (0), the 
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partial decay widths of different channels change in the same fashion. But the 
branching ratios remain unchanged. The striking feature is that, for some radion 
masses the gluon gluon branching ratio is almost equal to 1 . The other thing to 
note is that the photon photon branching ratio is also 5 to 8 times larger than 
the 2 photon branching ratio of the SM higgs. We see in the next section that 
this will have some interesting consequence in radion search. In this plot the SM 
higgs mass of 120 GeV is used for the purpose of illustration. This value of the 
SM higgs is well above the current experimental lower bound from the LEP Q . 
Another interesting feature of this plot is that the radion branching ratio into 
gg or 77 increases slowly with after the sharp fall around WW threshold. 
On the contrary the higgs branching ratio in these two modes decreases with 
nifi after the WW threshold. 



2 Radion Production at Hadron Colliders 

The radion coupling to weak gauge bosons are suppressed relative to hWW and 
hZZ couplings by a factor of j^. It is known that higgs boson production at 
LHC by the weak boson fusion mechanism is itself suppressed relative to the 
gluon fusion mechanism over most of the range of Mh {Mh < 1 TeV). Hence 
the radion production at LHC by weak boson fusion is also not expected to be 
the dominant or efficient mechanism. Also the radion couplings to light valence 
quarks are extremely small. Therefore in this article we shall focus on radion 
production by gluon fusion. The radion production cross-section at a hadron 
collider via gluon fusion mechanism is given by, 

a{p p{p) -.4> + X) = T{4>^gg)^r r — g{x) g{-) (7) 

<t> 

Here r is a dimensionless variable given by where S is the proton - 
proton (anti-proton) center of mass energy. Note that L^gg generates momen- 
tum dependent (j)gg, <j)ggg and 4>ggg couplings. The strength of these couplings 
are proportional to P{gs) which includes the contribution of gluons as well as 
the dynamical quarks. The Lagrangian also indicates that radions could be 
produced either singly or in association with gluons at a hadron collider by 
gluon fusion. In fact the gluon fusion mechanism turns out to be the dominant 
production process for radions at hadron colliders over most of the interesting 
range, qg ^ <j) and qg ^ (p does make a small contribution to radion produc- 
tion in 0(asY ■ The radion coupling to gluons is very similar in structure to 
the effective Lagrangian that gives the higgs coupling to gluons [5] in the heavy 
quark limit. As in the case of the radion the gluon fusion also turns out to 
be the primary production mechanism for higgs boson at hadron collider. The 
dominant contribution to gg — > h arises from closed loops of heavy quarks that 
occur in the theory. In this work we shall assume that the number of heavy 
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quarks (Nh) is equal to one namely the top. It has been shown |^ that the 
heavy quark limit ^ — > is an excellent approximation to the exact two loop 
corrected rate ioi gg ^ h. As mh ^ 2 rrit the exact result rises above the heavy 
quark result and exhibits a small bump corresponding to the ti threshold. The 
width of the bump i.e the departure region increases with increasing nit- The 
disagreement between the two results in the nih > 2mt region however is always 
less than a factor of two at LHC. The heavy quark limit for gg h can be 
obtained from the gauge invariant effective Lagrangian |Q 



1 _ 2 (3h 



G"'''' - — hit. (8) 



gsil + 5) V 

6 = 1 + 2^ is the anomalous dimension of the mass operator arising from 
QCD interactions. Ph is the heavy quark contribution to the QCD beta function. 
Since the hgg coupling in the Mq — + oo limit arises from heavy quark loops it 
is only the heavy quarks that contribute to the /3h in eqn(3). To order (a^) the 
heavy quark contribution [6] to (3{gs) is given by f3h = + 
the other hand the P{gs) that appears in the gluon coupling to the dilaton like 
radion mode arises from the trace anomaly. The trace anomaly has its origin 
in the heavy regulator fields as their masses are taken to infinity. So it includes 
the gluonic contribution as well as that of dynamical quarks. This difference in 
the two beta function contributions makes the (j)gg coupling greater than the 
hgg coupling even for = j. However with increasing (0) the hgg coupling 

ultimately wins over the (i)gg coupling. This feature is clear from fig. |^ where it 
is shown that with increasing (0), a{pp — s- (j)) ultimately becomes smaller than 
a{pp-^ h). 

Let us now make some rough numerical estimate about the ratio '^J'^^p^jh) 
in the lowest order {0{as)). Our estimates will depend only on the relative 
strength of 4>gg and hgg couplings. For \fi < 2mt, — — 3.84jj whereas 

/3ft, = to lowest order. Also in this region the heavy quark limit provides 
a good approximation to the exact result for aipp h). We find that « 

-2.88 for (0) = 1 TeV. Above the 2mt threshold we have « -3-50^. 

Although in this region the heavy quark limit does not work tfiat well for the 
higgs cross section we can still get an order of magnitude estimate (lower by at 
most a factor of two) using it. The ratio of couplings now > 2mt) becomes 
^29* w -2.63. 

9ggh 

Using the fact that the effective Lagrangian for higgs and radion production 
by gluon fusion are similar except for couplings we find that if ((/')=1 TeV then 
the cross section for pp ^ cf) will exceed that of — + ft- by a factor of 8.3 for 
y/§ < 2mt and by a factor of 6.9 for ^/l > 2mt. However for {4>)—5 TeV the 
radion production cross section will be suppressed relative to the higgs cross 
section roughly by a factor of three. These features have been exhibited in fig. 
H where we have plotted the lowest order higgs production cross section (both 
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Figure 2: Radion production cross-section (a,b,e) at the LHC, for different val- 
ues of {(j)) . We also present the SM higgs production cross-section (c) the exact 
result and (d) in the heavy quark limit. 



exact and the heavy quark hmit) and the radion cross section (for three different 
values of (0)) against the mass of the particle ((/> or h) at LHC. 

The estimates given above are based on lowest order calculations. It is 
known that higher order QCD corrections increases the lowest order rate by a 
factor (K factor) that lies between 2 and 3 at LHC ^. The QCD radiative 
corrections done in the heavy quark limit forms an excellent approximation to 
the exact calculations. To calculate the K factor one therefore always uses the 
heavy quark limit. But in the heavy quark limit the effective Lagrangian for 
higgs production is similar in structure to the Lagrangian for radion production. 
Hence the K factor for higgs production in the heavy quark limit will be the 
same for the radion also. So higher order QCD corrections will not affect the 
relative rate between the radion and the higgs to a very high degree of accuracy. 
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3 Detection and Possible SM Backgrounds 



Let us now concentrate on the detection of radion at a hadron collider like 
Tevatron. The dominant decay mode of a 50-150 GeV radion as can be seen from 
fig. 1^, is to bb or to gg. But the striking feature is that the 77 branching ratio of 
the radion is larger than the higgs case by a factor 5-8, over a considerable mass 
range. The higgs production rate by gluon fusion and and its decay into the 77 
mode are both suppressed relative to that of the radion. This is the reason why 
two photon final state is not a good bet for the higgs at the Tevatron. 

At the Tevatron radion production cross-section varies from 140 pb to 1 pb as 
we vary radion mass from 20 (current lower bound on radion mass comes from 
LEP -II §) to 160 GeV. We have included a NLO QCD correction factor of 
2 in our calculation. This cross-section with the presently collected luminosity 
will give rise to some 10 radions. If the radion decays into the 77 mode 
the final state will consist of 2 hard photons. The main background of this 
2-photon final state comes from the pair annihilation of the valence quarks 
and anti-quarks. The other dominant source of 77 background is gluon gluon 
annihilation to two photons. Though this is suppressed to the former by a factor 
of aj. , dominance of gluon flux over the quarks flux, can make this comparable 
with the former. We do not calculate this second contribution explicitly. We 
multiply the qq — > 77 contribution by a factor of 2 to take this into account. At 
the Tevatron this is a conservative approximation. We have used a parton level 
monte-carlo event generator to estimate the numbers for both the signal and 
the background and the CTEQ-4M parametrisation for the parton densities 
in our entire analysis. 

The following cuts have been applied to differentiate between signal and 
background. 
pI > 10 GeV 

We demand the photons are in the central part of the detector. 
1^1 < 3.. 

We also require that the angular separation between the photons be substantial, 
i.e. ARjj = ^(Ar/^^)2 + (A(/>^^)2 > 0.5 

Even after applying these cuts, the SM background cannot be removed com- 
pletely. So the strategy is to compare the invariant mass (of the photon pair) 
distribution of the signal and background. For the signal (of a definite m^), in- 
variant mass distribution shows a sharp peak over the continuum background. 
The sharpness of the peak depends mainly on the detector resolution and the 
partial width of ^ 77. In fig. ^ we plot the invariant mass (M^^) distribution 
for the signal (dots) and backgrounds (broken histogram) assuming an uniform 
bin size of 5 GeV. Here one can see that, for low M^^, number of background 
events is higher than the signal. But as M^^ increases, the number of back- 
ground events in the mass bins falls off more rapidly than the signal- which 
more or less remains the same over the entire mass range we are interested here. 
Once the radion mass is near to 2mw, signal events falls off sharply, due to the 
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Figure 3: Invariant mass distribution for signal and background at the Tevatron 
for {(j)) = 1 TeV. 



sharp fall of 77 branching ratio. 

We find that the 77 mode is good enough to exclude radion mass nearly 
upto 120 GeV even with the presently collected luminosity of 110 pb~^ . Here 
we have taken the radion vev to be 1 TeV. On the other hand at the Tevatron 
with the presently collected luminosity one cannot say anything about the SM 
higgs. 

Now let us go over to the case of Upgraded Tevatron with center of mass 
energy of 2 TeV and luminosity of 1 fb~^. As center of mass energy is almost 
equal to the present, the signal and background are just 10 times larger than 
the previous case. This is evident from the fig. 4. In this figure we also plot 
the significance (= No.ofSignal \ signal upto M^j equal to 100 

Y No. of Background 

GeV. One can check, over this entire mass range, significance is greater than 

5. Generally a significance greater than 5 points to the discovery. One can also 
easily check that once M^-y is greater than 100 GeV, no. of background events 
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in the corresponding bins become less than 1 events. So in this region we can 
demand 5 signal events as a benchmark for discovery. Thus at the upgraded 
Tevatron one can discover radion mass upto 160 GcV and exclude upto 165 
GeV. Once M^^ is greater than 165 GeV, our signal falls off very sharply and 
we cannot say anything more about it. 
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Figure 4: Invariant mass distribution for signal and background at the Tevatron 
Upgrade for {<p) = 1 TeV. 

If we change {(p) the branching ratio of the radion to different channels 
remains same. So if we take {(j)) = 2 TeV, the cross-section and number of 77 
events become ^ of the present case ((0) = 1 TeV). And at the Tevatron with 
the presently collected luminosity, we cannot say anything about it. At the 
Tevatron Upgrade, we also cannot talk about the discovery, for <45 GeV. 
At higher masses, one can discover upto 150 GeV if (</>) = 2 TeV instead of 
1 TeV. 

Finally we want to examine the search prospects of this scalar particle at the 
LHC. Though the two photon branching ratio drops very sharply near ~ 
140 GeV, in contrast to the SM case it remains constant with mi after 140 
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GeV. Also the radion production rate is almost 8 times higher than that of 
higgs production rate. So unlike the SM higgs boson, for > 140 GeV, 77 
mode is a viable avenue to discover or exclude the radion for masses well upto 
1 TeV. 
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Figure 5: Invariant mass distribution for signal and background at the LHC for 
(0) ^ 1 TeV. 

Conservatively, we take the luminosity to be equal to 30 fb^^. The source of 
SM backgrounds remain the same here. But unlike the Tevatron, the 2-photon 
background coming from pair annihilation of quark and anti-quark becomes 
less severe. This is because the anti-quark coming from one proton has to be 
excited from the sea. But gluon density in the proton is much larger than the 
quark density at LHC. So to take into account the gluon gluon contribution to 
2-photon background we multiply the quark-anti-quark contribution by a factor 
of 6 [[ The cuts we used here are rather similar to the Tevatron case. Only we 

^ Generally at the LHC energies the gluon-gluon contribution to two photon background is 
5 times larger than the quark anti-quark contribution. For ruj-y = 110 (130) GcV this factor 
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change the cuts on the transverse momenta of the photons to be greater than 
20 GeV. We present in fig. 5 the number of signal events and the corresponding 
number for the background against the 77 invariant mass. To plot the invariant 
mass distribution for the background we choose a bin size of 5 GeV. Form the 
figure one can easily see that the background is an order of magnitude smaller 
than the signal over the entire mass range. And once M^^ is greater than 250 
GeV, there are less than one event in the respective mass bins. That's why we 
do not show it in the figure. So LHC can easily discover radions with masses 
upto 1 TeV for (0) = 1 Tev. If we increase {(j)) to 4 TeV one can easily check 
from fig. 5 that our discovery limit comes down to 650 GeV. 



4 Conclusions 

In this paper we examined the radion production and its subsequent decay into 
SM particles at hadron colliders. The radion production cross-section is larger 
than the higgs production cross-section by a factor of 6-8. The partial width of 
the radion to two gluons and two photons is also enhanced due to the enhanced 
(f> g g or (j) -f J coupling. The other partial widths are suppressed with respect 
to higgs widths by a factor that depends on {(j>) . We also discussed the viability 
of the 2-photon signal for the radion at the Tevatron. One can exclude radion 
mass upto 120 GeV even with the presently collected luminosity. Upgraded 
Tevatron with higher center of mass energy and higher luminosity certainly can 
discover radion upto 160 GeV mass, if not it can exclude it upto 165 GeV. 
Similarly at the LHG one can definitely discover radions with masses up to 1 
TeV. These estimates are done assuming (0) =1 TeV. With increasing (ip) both 
the exclusion limit and the discovery limit using the 77 mode will come down. 



Note added: While this work was being completed there appeared 
one paper [12| which discusses some of the issues presented here. 
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